Fermi surface renormalization and confinement in two coupled metallic chains 
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Using a non-perturbative functional renormalization group approach involving both fermionic and 
bosonic fields we calculate the interaction-induced change of the Fermi surface of spinless fermions 
moving on two chains connected by weak interchain hopping t±. For a model containing interband 
backward scattering only we show that the distance A between the Fermi momenta associated with 
the bonding and the antibonding band can be strongly reduced, corresponding to a large reduction 
of the effective interchain hopping oc A. A self-consistent one-loop approximation neglecting 
marginal vertex corrections and wave-function renormalizations predicts a confinement transition 
for sufficiently large interchain backscattering, where the renormalized vanishes. However, a 
more accurate calculation taking vertex corrections and wave-function renormalizations into account 
predicts only weak confinement in the sense that < l^xl- Our method can be applied to 

other strong-coupling problems where the dominant scattering channel is known. 

PACS numbers: 71.10.Pm, 71.27.-f a,71.10.Hf 



I. INTRODUCTION 

In strongly correlated Fermi systems electron-electron 
interactions can have drastic effects on the geometry and 
the topology of the Fermi surface. For example, strong 
forward scattering can give rise to a Pomeranchuk insta- 
bility, where the shape of the Fermi surface spontaneously 
changes such that it has a lower symmetry than the un- 
derlying lattice"^. Another example is the Lifshitz tran- 
sitioni, where the topology of the Fermi surface changes 
discontinuously without symmetry breaking as a func- 
tion of some external control parameter. This gives rise 
to anomalies in thermodynamic and kinetic properties 
of a metal. Conditions on the range and the strength 
of the interaction leading to Pomeranchuk and Lifshitz 
transitions have recently been derived in Ref . [3] . 

In this work we shall focus on another type of phase 
transition associated with the geometry of the Fermi sur- 
face, which we call confinement transition. This quantum 
phase transition can occur in quasi one-dimensional met- 
als with an open Fermi surface, consisting of two discon- 
nected weakly curved sheets. Due to strong interactions, 
the curvature of the Fermi sheets can be smoothed out 
and can eventually vanish in certain sectors. In the ex- 
treme case, the renormalized Fermi surface consists of 
two completely flat parallel planes. The motion of the 
fermions in real space is then strictly one-dimensional, 
although in the absence of interactions it is not. We 
therefore call such a transition confinement transition. In 
the confined state the low-energy properties of the sys- 
tem resemble that of a one-dimensional Luttinger liquid. 
Because the Fermi surface in the confined state has an 
additional nesting symmetry, at the confinement transi- 
tion the symmetry of the Fermi surface increases, in con- 
trast to the Pomeranchuk instability, where interactions 
lower the symmetry of the Fermi surface. An interaction- 
induced flattening of the Fermi surface might also play 
a role in the Hubbard model close to half flUing, where 



the bare Fermi surface consists of four almost flat sec- 
tors. Completely flat parts of the Fermi surface can give 
rise to non-Fermi liquid behaviori^i^. Evidence of an 
interaction-induced flattening of the Fermi surface of the 
Hubbard model close to half flUing has been found in 
Ref. a. 

Similar to the Pomeranchuk transition, the conflne- 
ment transition is a strong-coupling phenomenon. Hence, 
the usual weak coupling perturbative expansions are 
not sufficient to study this transition. Due to a lack 
of controlled methods to deal with strongly interacting 
fermions in dimensions larger than one, it is very diffl- 
cult to study the confinement transition. To shed some 
light on the underlying mechanism, we shall in this work 
consider the simpler problem of just two metallic spin- 
less chains coupled by weak interchain hopping tj_. The 
confined state corresponds to a vanishing renormalized 
interchain hopping t*j_ = 0, so that electrons cannot move 
from one chain to the other, in spite of the fact that the 
bare interchain hopping is finite. In a subsequent paper—, 
we shall discuss the more difficult confinement problem in 
an infinite array of coupled chains. It turns out that the 
basic mechanism responsible for the tendency towards 
confinement can already be understood from the simpler 
two-chain problem. 

Because perturbation theory in the two-chain prob- 
lem is plagued by the usual infrared divergencies of one- 
dimensional Fermi systems, even in the limit of weak in- 
teractions the Fermi surface cannot be calculated within 
renormalized perturbation theorjs^iiiS. In dimensions 
D > 2 the Fermi surface deformation has been studied 
to all orders in perturbation theory in Ref. [ll|| for a gen- 
eral class of models. Within the framework of the renor- 
malization group (RG) the Fermi surface can be defined 
non-perturbatively from the requirement that the rele- 
vant coupling constants ri{kp) related to the self-energy 
T,{kp,uj = 0) at the true Fermi surface kp flow into a 
fixed point of the R C^^i^^ . In Ref. we have calculated 



2 



the shift of the Fermi surface in the two-chain system 
within the usual weak coupHng expansion of the RG (3- 
functions. We have shown that interchain backscattering 
gives rise to the dominant logarithmic renormalization of 
the distance A = k^ — k~ between the Fermi momenta fc+ 
and associated with the bonding and the antibond- 
ing band. Denoting by Ai the value of A within the 
Hartree-Fock approximation, the self-consistency condi- 
tion for the true Fermi point distance in the spinless two- 
chain system can be cast into the form 



A 



Ai 



1 + 2.g^ ln(Ao/A) 



(1.1) 



where Aq is an ultraviolet cutoff and go is the bare value 
of the dimensionless coupling constant describing inter- 
chain backscattering, which will be defined more precisely 
in Sec. lIVBl From Eq. (jl.ip we see that sufficiently large 
interchain backscattering strongly reduces the value of A. 
But the renormalized A never vanishes, so that there is 
no true confinement transition. One should keep in mind, 
however, that Eq. has been derived by expanding 

the RG /3-functions to second order in the coupling con- 
stants, so that it is not allowed to extrapolate this ex- 
pression to large values of 170 ■ 

To find out whether in the spinless two-chain system 
sufficiently strong interchain-backscattering can give rise 
to a confinement transition where the renormalized ef- 
fective interchain hopping oc fc"*" — k~ vanishes, we 
use here a generalization of the collective field functional 
RG approach with momentum transfer cutoff developed 
in Ref. (isj . It turns out that with this approach we can 
analyze the regime where the dimensionless interchain 
backscattering interaction is of the order of unity. The 
crucial point is that from the weak coupling analysis^^ we 
know that the confinement transition is driven by inter- 
chain backscattering, so that it is natural to decouple 
the interaction in this scattering channel via a suitable 
bosonic Hubbard-Stratonovich field. Simple truncations 
in the resulting mixed boson-fermion theory correspond 
to infinite resummations in an expansion of the RG /3- 
functions in powers of go. 



II. EFFECTIVE LOW-ENERGY MODEL 

We consider spinless fermions confined to two chains 
that are coupled by weak interchain hopping t±. The 
kinetic energy of the two-chain system is diagonalized by 
forming symmetric (bonding band) and antisymmetric 
(antibonding band) combinations of the eigenstates as- 
sociated with isolated chains. Denoting by the energy 
dispersion of a single chain in the absence of interchain 
hopping, the energy dispersion of the non-interacting sys- 
tem is = efc — at±, where a = +1 labels the bonding 
band and cr = — 1 labels the antibonding band. It is use- 
ful to think of cr as a pseudospin labcl^^<, in which case 
t± = h corresponds to a uniform magnetic field h in the 
z-direction. 



The problem of finding the low-energy properties of 
two coupled metallic chains has been studied previously 
by many author j^"-^'^-^^-^'''--^^-^^-^"-^^-^^'^^-^'^-^^-^^'^'^-^^-^^ 
However, the problem of self-consistently calculat- 
ing the true Fermi surface has only recently been 
addressediSiiiiiS^. At low energies all possible scattering 
processes in the spinless two-chain system can be divided 
into four different classesi^: (1) forward scattering pro- 
cesses, parameterized in terms of three different coupling 
constants f~^^, f and /"' — f '", where the labels 
indicate the band of the fermions involved in the scat- 
tering process; (2) interband backward scattering, which 
in pseudospin language corresponds to transverse spin- 
exchange, so that we shall call the corresponding dimen- 
sionful coupling constant J"*" (the associated dimension- 
less coupling go will be introduced in Sec lIVB]) : (3) pair- 
tunneling, which can also be viewed as interband Umk- 
lapp scattering, parameterized in terms of a coupling con- 
stant by u; and finally (4) intraband Umklapp scattering, 
which is expected to be important only at commensurate 
fillings. Neglecting the latter process and setting for sim- 
plicity / , the low-energy interactions can be 
expressed in terms of four marginal coupling constants 

/ = |(/+- + /++), ^" = - /++), J^, and u. 

In the bonding-antibonding basis the system can then 
be modeled by the following effective Euclidean action in 
pseudospin notation. 



1 r r 



K 



[fik)pKPK ~ jHk)mKmji 



Jk 

(2.1) 

where = ej, — /i — ah, and we have introduced the 
following composite fields. 



K+K 



JK 



K+K 



^K 



K 



'I'K^K+K 



(2.2a) 
(2.2b) 
(2.2c) 



We use the imaginary time formalism at zero tempera- 
ture and have introduced collective labels K = {k, iuj) for 
fermionic fields and K = {k, iuj) for bosonic fields, with 
the notation Jj^ = / j^^- Note that the Fourier com- 
ponents of the density and the longitudinal spin-density 
field satisfy p-K = Pk and tti-k — rnx, while the spin- 
fiip field sk is complex and do not have this symmetry. 
The interaction functions f{k), J"(A:), J'^(k), and uik) 
should be considered as phenomenological low-energy 
couplings which are only non-zero for |fc| < Aq <C kp. 
Hence, these couplings should not be directly compared 
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with the bare couphng constant in the Hubbard model^S. 
The signs and normahzations in Eq. (|2.ip are chosen 
such that for = the model has rotational invari- 
ance in pseudo-spin space, and that for the Hubbard 
model all couplings are positive^*^. However, in our ef- 
fective low energy model there is no reason to expect 
rotational invariance in pseudospin-space, so that in gen- 
eral Jll(fc) 7^ J-L(fc). 

The model defined in (|2.1|) is still quite complicated 
and contains many interaction processes which are not 
essential for the confinement transition. In fact, from our 
previous weak coupling analysis^* we know that the dom- 
inant renormalization of the difference between the Fermi 
points is due to the interchain backscattering process de- 
scribed by the coupling J^(fc). In this work we study 
a minimal model describing the confinement transition 
by simply neglecting the forward scattering interactions 
fik) and J"(fc), as well as the pair tunneling coupling 
uik) in Eq. (|2.ip . However it is known^^'^^ that suffi- 
ciently strong pair tunneling can destabilize the Luttinger 
liquid phase; we shall come back to this point in Sec. |Vl 
In pseudospin language, our model then describes a one- 
dimensional spin S* = 1/2 Fermi system subject to a uni- 
form magnetic field in z-direction with an attractive fer- 
romagnetic spin exchange involving only the transverse 
{xy) spin components. The latter tends to align the spins 
in the direction perpendicular to the magnetic field. The 
Fermi surface renormalization is essentially determined 
by the competition between the a;?;-exchange interaction 
and the external constraint imposed by the uniform mag- 
netic field, which tends to align the spins along the z-axis. 
The phase diagram of the model (|2.ip in the space of all 
couplings has been discussed by Fabrizic^^. The qualita- 
tive behavior of the weak coupling RG flow in the space 
of couplings spanned by j" , J-^ and u is shown in Fig. [T] 
Obviously, there is a finite regime in coupling space where 
the spinless two-chain system is a stable Luttinger liquid, 
with gapless excitations. In this work we shall assume 
that the qualitative fixed point structure suggested by 
the weak coupling analysis remains correct even in the 
strong coupling regime. We can therefore choose the bare 
parameters such that the system belongs to the basin of 
attraction of the Luttinger liquid fixed point manifold. 

At low energies we may further simplify our model (at 
least in the deconfined phase) by linearizing the energy 
dispersion at the Fermi surface, which for our model con- 
sists of four points ak'^ , where the chirality index a = ±1 
labels the left/right Fermi point. Note that the true 
Fermi points are defined via 

eak' - ^J^-(yh + T.''{ak'',i{))^Q , (2.3) 

where I]'^(aA;°', iO) is the exact self-energy for vanishing 
frequency and for momenta at the true Fermi surface 
ak'^ of the interacting system. To linearize the energy 
dispersion at the true Fermi surface, we add and subtract 
from the non-interacting energy dispersion the counter- 
term 

^ -Y.^iak" ,iQ) , (2.4) 




FIG. 1: (Color Online) Qualitative behavior of the weak cou- 
pling RG flow of the model (|2.1|l in the space of coupling 
constants j", and u. The thick black line is the line of 
fixed points describing the stable Luttinger liquid phase. 



and approximate 

= eak-+k - fi-ah + Y.^iak", iO) - Y" {ak'^ , iO) 
^av^k + fi^, (2.5) 

where Vq is the bare Fermi velocity at the true Fermi 
surface for pseudospin a. In analogy with the definition 
of the couplings in the Tomonaga-Luttinger model^, we 
now generalize the interaction by introducing chirality in- 
dices — > J^^, . We shall refer to the diagonal processes 
J^Q, as chiral interactions (these are called 54 processes in 
the Tomonaga-Luttinger model). Similarly, off-diagonal 
elements J a -a be called non-chiral processes (corre- 
sponding to the (72-pi'Ocesses in the Tomonaga-Luttinger 
model). Defining new fields 

we replace the action (|2.ip by the following effective low- 
energy action describing the physics of the confinement 
transition in our system of two spinless metallic chains, 

aa' 

where it is understood that the /c-integration has a mo- 
mentum transfer cutoff |fc| < Aq <C kp, and 

Si?a = J^^Kc^Uko. ■ (2-8) 
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FIG. 2: Bare interaction vertices of the action Si[^] given 
in Eq. (|3.4|l . The fermionic fields ip'^ and 1/)°^ are denoted by 
sohd arrows, with the spin-projection a = ±f written next 
to the arrows. Bosonic spin-flip fields x a^nd x are denoted 
by wavy arrows. Incoming arrows denote 1/''^ and X) while 
outgoing arrows correspond to the conjugate fields tj}'^ and X- 

III. EXACT RG FLOW EQUATIONS 

A. Hubbard Stratonovich transformation 

Because the confinement transition is a strong cou- 
pling phenomenon, the usual weak coupling RG approach 
based on the expansion in powers of J^^, is not suitable. 
To develop a RG approach which does not rely on a weak 
coupling expansion, we decouple the spin-flip interaction 
with the help of a complex Hubbard-Stratonovich field 
Xa- For convenience we collect all fields in a super-field, 

•J" = V'a, V'a,X«,X«) • (3-1) 

Taking into account that there are two chiralities a = ±1, 
our super-field has totally 12 components (8 fermionic 
and 4 bosonic ones). The ratio of the partition functions 
with and without interactions can then be written as 

Z _ /P[$]e-5o[*l-Si[*] 

Z'o " /I?[$]e-So[*] ' (^-^^ 

with the Gaussian part of the effective action given by 

+ jA-^XlxKcXRa' • (3.3) 

Here is a matrix in chirality space with matrix ele- 
ments given by J^^, , and the interaction is 

= E / [^KaXKa+ Sji^XRo] ■ (3-4) 

A graphical representation of the bare interaction vertices 
in Eq. ((^ is shown in Fig. [H The coupled RG flow 
equations for the one-line irreducible vertices of the above 
mixed boson-fermion theory can be obtained as a sp ecial 
case of the general flow equations given in Ref. [15l |. 

B. Functional RG flow equations in momentum 
transfer cutoff scheme 

In order to calculate the true Fermi surface, we need to 
know the exact counter-term /Iq = — S(afc°',iO), which 



can be obtained from the flowing self-energy E^(ii', a) 
in the limit of vanishing infrared cutoff A — > 0. The 
form of the flow equation for S^(iir, a) depends on 
the RG method employed. Here we use the hierarchy 
of functional RG equations for the one-line irreducible 
vertices22ii^ of mixed boson-fermion models developed 
in Ref. [T^ . A similar approach involving both fermionic 
and bosonic fields has been developed in Refs. 34,3^. 
In principle, one can also obtain the flowing self-energy 
within the purely fermionic parameterization of the hier- 
archy of flow equation9i2iiiii2&. However, with the usual 
truncations necessary in this approach it is not possible 
to reach the strong coupling regime. 

In the momentum transfer cutoff schemeJ^ we impose 
a cutoff A only on the momentum k transfered by the col- 
lective bosonic field. The resulting RG flow equation for 
the fermionic self-energy (K, a) is shown graphically 
in Fig. [21 The corresponding analytic expression is 

= [ Ff{K,a)Tf^\Ka,-Ka;K,-K,a) 

J K 

+ / Ff{K,a)Gl{K + K + aa 1:^,0) 

JK 

xF^^'^^ (Kcr; K + K,a; -K, a) 

xT^^'^\k + K,a;K,a;K,a) . (3.5) 

Here {K, a) is the flowing fermionic single-particle 
Green function for a given pseudospin a and chirality 
index a. We use the notation a — —a and measure the 
wave-vectors k with respect to the true Fermi momenta 
afc'^, defining 

Gl{K, a) = Gliak" + k, ito) , (3.6) 

and 

Gl{K + aaA, a) = Gl{ak'^ + k + aaA, iuj) . (3.7) 

The shift aaA = aa{k^ ~ k~) in the argument of in 
Eq. p.Sp is due to the fact that in G^ (K, a) the wave- 
vector k is measured with respect to a different refer- 
ence point than in G']^{K,a). The function F^°'{Ka) in 
Eq. (|3.5p is the single scale bosonic spin-flip propagator, 
which is defined by 

FfiKa) = ^S{\k\ - A)[Fr{K)U , (3.8) 

where F^*^ (K) is a matrix in chirality space whose inverse 
has the matrix elements 

[Fr{K)Va' - [2^X1 - Saa'HriK, a) , (3.9) 

where Il'J^°'{K,a) is the flowing spin-flip susceptibility. 
In the momentum transfer cutoff scheme, the RG flow 
of H^°"(^,a) is driven by the one- line irreducible vertex 
with four external boson legs, as shown in Fig. 21 The 

vertices r^^'^\Ka;K'a;K,a) in Eq. jHS]) are the flow- 
ing spin-flip vertices with two fermionic and one bosonic 





FIG. 3: Exact flow equation for the fermionic self-energy 
EX (AT, a) in the momentum transfer cutoff scheme. The thick 
sohd arrow is the flowing fermionic Green function and the 
thick wavy hue with a slash is the flowing single scale spin- 
flip propagator deflned in Eq. (|3.8|) . The one- line irreducible 
vertices are represented by shaded triangles. A label (n, m) 
inside a shaded triangle means that the vertex has n fermionic 
and m bosonic external legs. 



external legs. In the momentum transfer cutoff scheme 
these vertices satisfy the exact flow equations shown in 
graphically in Fig. [5l with initial condition 



(3.10) 



(2 2) 

Finally, the vertex T)^ on the right-hand side of 
Eq. (|3.5p is the one-line irreducible vertex with two 
fermionic and two bosonic external legs. We do not give 
the flow equation for this vertex, because purely bosonic 
vertices with more than two external legs and mixed ver- 
tices with two fermionic and more than one bosonic ex- 
ternal leg have negative scaling dimensions and are ir- 
relevanti^. We expect that their effect can be implicitly 
taken into account by re-defining the numerical values of 
the relevant and marginal couplings 

The initial condition for the fermionic self-energy at 





FIG. 4: Exact flow equation for the spin-flip susceptibility in 
the momentum transfer cutoff scheme. 



scale A = Ao is simply 



i]X„(i^,«) = o 



(3.11) 



Similarly, the vertices with two fcrmion legs and more 
than one boson leg appearing on the right-hand side 
of the flow equation for the spin-flip vertex shown in 
Fig. [5] also vanish at the initial scale. However, the 
price we pay for introducing a cutoff only in the mo- 
mentum transfer are non-trivial initial conditions for the 
purely bosonic vertices, which are initially given by closed 
fermion loopsi^. In particular, the loop with two exter- 
nal boson legs is initially given by the non-interacting 
spin-flip susceptibility. 



K 



Gl^iK,a)Gl^{K + K + aaA.a) 



(3.12) 



where for our model with linear energy dispersion, 

1 



ibj — avQk — /Xg 



Denoting by 



Ao 



(3.13) 



(3.14) 



the distance between the Fermi momenta fcg and in 
the absence of interactions, the relation between the true 
distance A = k'^ — k~ and Aq can be expressed in terms 
of the counter-terms /ig — — S(q;A:°', iO, a) as follows, 



(3.15) 



see also Eq. (p:iT|) below. Usi ng Eq s. and ([XTC]) 

we can explicitly evaluate Eq. (|3.12[) . 



(crAo + ak) 



2ttvq VQ{aAo + ak) — ioj 



(3.16) 



C. Rescaled flow equations and classification of 
vertices 

To classify the various vertices according to their rele- 
vance, it is useful to make them dimensionless by multi- 
plying them with a suitable power of the running cutoff 
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2,3 




+ 
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+ 




FIG. 5: Exact flow equation for the spin-flip vertices in the 
momentum transfer cutoff scheme. 



A. Following Ref. [l^, we define dimensionless fermionic 
labels Q — {q,ie) = {k/A,iuj/n\), and bosonic ones 
Q = {q,H) ~ {k/A,iuj/D,A_). Here fl^ — vpA, where vp 
is some average Fermi velocity. For simplicity we shall 
write the above relations as Q = K/A and Q = K/A. 
We consider all rescaled quantities as functions of the 
logarithmic flow parameter I — ln(Ao/A). 

In order to define the Fermi surface within the 
framework of the RG, we subtract the counter-term 
I]°'{ak'^ ,iO,a) = — from the flowing self-energy and 



then rescalei^, 



Z? 



SZa 

= §^[EX(AQ,a) + Ai[J] . 

"A 



(3.17) 



Here is the flowing wave-function renormalization fac- 
tor, which is defined in terms of the flowing self-energy 
as follows, 



(3.18) 



E=0 



The corresponding rescaled fermionic propagator is 

Gr(Q,«)-^GX(AQ,a). (3.19) 

The rescaled self-energy satisfies the exact RG flow 
equatio n^^i^^ 

9,Sr(Q, a) = (1 - ifi - qd, - ed,)t^{Q, a) + f ^(Q, ") , 

(3.20) 

where the flowing anomalous dimension of the Fermi 
fields is 



-di\nZ[ 



dtf{0,ie,a) 



d{te) 



(3.21) 



e=0 



and the function Tf{Q, a) follows from Eq. p. 51 
tl{Q,a) 



^[-AdA^l{K,a) 



PriQ, a)f p'^' (Q, a; -Qa; Q, -Q, a) 



FriQ, a)GfiQ + Q + aaA^a) 



xr[^'^\Q,a-Q + Q,ij;-Q,a) 



where 



A _ A 

A "a;;' 



(3.22) 



(3.23) 



is the rescaled true difference between the Fermi points. 
The rescaled bosonic single scale propagator follows from 
Eq. 

Fr{Q,a) = -^AFr(AQ,a) 

^ S{\q\-l)[Fr{Q)U, (3.24) 

with 

[FriQ)Vo.' = ^/[2i^oJ^]-^, -<5a.'Iir(Q,«) , (3.25) 
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and the rescaled spin-flip susceptibility 



(3.26) 



Here Zi is the wave-function renormaUzation factor as- 
sociated with the bosonic spin-flip field Xcn and the con- 
stant with units of inverse velocity has been introduced 
to make all rescaled vertices dimensionless. The rescaled 
spin-flip vertex in Eq. p.22p is 

ff-^\Q,a-Q',a;Q,a) = [Zi/v^f^ [KZ+ /n^l^'^ 
xrf (AQ,(T;AQ',a;AQ,a) , (3.27) 

and the rescaled vertex with two external fermion and 
two boson legs is 

xrf (AQ',(T;AQ,CT;AQ',Ag,a) . (3.28) 

For convenience we now choose = A/^a = l/^^F so 
that the prefactor on the right-hand side of Eq. (|3.27|) 
reduces to [ZiZ+Z^Y/'^. 

Let us now classify the various vertices according to 
their relevance. First of all, the key quantity to obtain 
the counter-terms /Xg is the momentum- and frequency 
independent part of the rescaled self-energy Sf (Q, a) de- 
fined in Eq. (|3.17p . which we call 



rr = Er(0,a) = ^[SX(0,a) 

"A 



The couplings r'^ satisfy the exact flow equation 

airr = (i-r/nrr+rr(o,«), 

with initial condition 



Mo 



S'"(aA:'^,iO) 
w_fAo 



(3.29) 



(3.30) 



(3.31) 



There are two marginal couplings related to the self- 
energy. The first is the wave-function renormalization 
factor Zf , which according to Eq. (|3.2ip is related to the 
flowing anomalous dimension via 



diZ^ = -vilZt . 



(3.32) 



The second is the dimensionless Fermi velocity renormal- 
ization factor 



i1 = 



dT.'[{q,iO,a) 



d{aq) 

which satisfies the exact flow equation 
dtl{q,iO,a 



(3.33) 



q=0 



d{aq) 



(3.34) 



If we retain only relevant and marginal couplings, the 
rescaled fermionic propagator with energy dispersion lin- 
earized at the true Fermi surface is given by 



Gi{q,ie,a) 



1 



le — avi q — r\ 



(3.35) 



Apart from Z^ and , the third marginal coupling of 
our model is the momentum- and frequency-independent 
part of the rescaled spin-flip vertex defined in Eq. (|3.27p , 



7, =ff^'^(0,a;0,(7;0,a) . 
It satisfles a flow equation of the form 



(3.36) 



(3.37) 



g=0 



where t]i — —di In Zi is the flowing anomalous dimen- 
sion of the spin-flip field, and the inhomogcneity F; ' 
depends on the irrelevant higher interaction vertices in- 
volving more than one external boson leg shown in Fig.O 
In particular, from the right-hand side of Eq. (|3.28p it is 
clear that the vertex F^ ' with two external fermion and 
two boson legs is irrelevant with scaling dimension — 1. 
This and the higher order irrelevant vertices vanish at 
the initial scale Ag and we shall set them equal to zero, 
expecting that their effect can be implicitly taken into 
account via a redefinition of the numerical values of the 
relevant and marginal couplings^. An exception is the 
vertex F^^'^'' involving four external bosonic legs, which 
according to Fig. [4] drives the flow of the spin-flip sus- 
ceptibility in the momentum transfer cutoff scheme. In 
contrast to the other irrelevant vertices, the vertex F^^'^'' 
is finite at the initial scale A = Ag, where it reduces to a 
symmetrized closed fermion loop^^. Below we shall pro- 
pose a simple approximate procedure to take the renor- 
malization of the spin-flip susceptibility generated by this 
vertex into account. Finally, we note that the vertex with 
four external fermionic legs is also marginal, but in the 
momentum transfer cutoff scheme it does not directly 
couple to the flow of the fermionic self-energy. 



D. Defining the Fermi surface within the 
functional RG 



The general method to obtain the counter-terms nec- 
essary to construct the true Fermi surface within the 
framework of the functional RG has been developed in 
Refs. [miliil. Let us briefly recall the main idea. As 
long as the flowing anomalous dimensions -qf of the Fermi 
fields remains smaller than unity for Z — > oo, we may de- 
fine the true Fermi surface self-consistently from the re- 
quirement that the relevant couplings associated with 
the fermionic self-energy approach finite limits for I ^ oo. 
This requires fine tuning of the initial values Tq , which 
defines a relation between rg and the flowing couplings 
on the entire RG trajectory. In higher dimensions, where 
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the Fermi surface is a continuum, infinitely many rele- 
vant couplings ri{kp) have to be fine tuned to define the 
Fermi surface. In the usual classification of critical fixed 
points, the Fermi surface thus corresponds to a multicrit- 
ical point of infinite order. Once the proper initial values 
Tq are known, the exact self-energy Yi" {ak'^ , iQ) can be 
constructed using Eq. (|3.3ip . 

Yr{ak'',iQ)^-^ll^-VF^ar'^■ (3.38) 

The requirement that flows into a RG fixed point im- 
plies for the initial valuesi^, 

Jo 

where 



(3.39) 



dtVt 



(3.40) 



is the average of the flowing anomalous dimension 
along the RG trajectory, and we have written (0) = 
rj^(0,iO,Q!) to emphasize that this quantity is actually 
independent of the chirality index a. For our effective 
model with linear energy dispersion we obtain for the 
Fermi point distance at constant chemical potential [see 
also Eq. (|231l ]. 



Ao 



'0 '0 



■^0 '"0 
noo 



poo 

^Ao^J2- dle-^'-'''^'tnO,a) , (3.41) 



where we have defined 



A = — , An - " ° 



Ao 



Ao 



(3.42) 



IV. 



CALCULATION OF THE TRUE FERMI 
SURFACE 

A. Truncation based on relevance 



Because a possible confinement transition is expected 
to be a strong-coupling phenomenon, the usual perturba- 
tive weak coupling RG^^ is not sufficient. We therefore 
propose an alternative truncation scheme based on the 
truncation according to relevance in the RG sense. In 
our model we have to keep track of the RG flow of the 
two relevant couplings rf , a = ±1, and the marginal 
couphngs Z[, vf , and ji. The couphngs rf , Z[ and 
associated with the fermionic Green function satisfy the 
flow equations given in Eqs. (|3.30|[3?32ll3.34p . The func- 
tion Ti{Q,a) appearing on the right-hand side of these 
equations is in general given in Eq. p.22p ; approximating 



the fermionic Green function by Eq. (|3.35p and the spin- 
flip vertex by its momentum- and frequency independent 
part 7/, we obtain 



ti{q,ie,a) 



dqde 
(2^ 



m\ - 1) 



i(e + e) - avf {q + q) - crijf A,* - rf 



(4.1) 



Here Af is the rescaled true difference between the Fermi 



points defined in Eq. (|3.23p . and the rescaled bosonic 
spin-flip propagator Ff"'{Q) is defined in Eq. (|3.25p . 

To calculate the Fermi surface, we need additional flow 
equations for the marginal part of the spin-flip vertex 
7/ and for the flowing spin-flip susceptibihty U1'^{Q,a). 
As far as 7; is concerned, we note from Fig. \5\ that the 
inhomogeneity r[^'^'' in Eq. (|3.37p which drives the flow 
of 7/ involves vertices with two fermionic and more than 
one bosonic external legs. These vertices are irrelevant 
and vanish at the initial scale Aq, so that it is reasonable 
to neglect them. We therefore set Fj — 0. We shall 
also neglect the bosonic wave-function renormalization, 
setting Zi = 1. In this approximation the flow of the 
rescaled spin-flip vertex is driven by the fermionic wave- 
function renormalization, 



dm 
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(4.2) 



Before discussing the spin-flip susceptibility Ilf^ (Q, a) , 
note that Eqs. (|3.34p and (|4.ip imply for the Fermi ve- 
locity renormalization factor 



which yields for the difference 



l^-ivt + Vi )(«/ 



(4.3) 



(4.4) 



Keeping in mind that r^f > 0, Eq. (|4.4p implies that a 
small initial difference between the Fermi velocities de- 
creases under renormalization. Thus, if the initial dif- 
ference — Vq is small and negligible, it becomes even 
smaller as we iterate the RG. Since the flow of the other 
couplings is not sensitive to a small difference in the 
Fermi velocities, it is consistent to approximate Vq « vp, 
so that from now on we shall set vf = 1. 

To close our system of flow equations, we need an equa- 
tion for the rescaled spin- flip susceptibility 11^'^ (Q, a), 
which in turn determines the flow of the spin-flip propa- 
gator as given in Eq. (|3.26p . In the momentum transfer 
cutoff scheme, the flow of flf^ {Q, a) is driven by the one- 
line irreducible vertex rj^^'^"* with four external bosonic 
legs, as shown in Fig. 21 Although this vertex is irrele- 
vant, it is flnite at the initial scale Ao, in contrast to the 
higher order vertices that drive the flow of the spin-flip 
vertex 7; shown in Fig. [5l It is therefore important to 

at 



take the renormalizations of 11^'^ (Q, a) due to F 



(0,4) 
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least approximately into account. Guided by the initial 
condition (|3.16p for the spin-flip susceptibility, we pro- 
pose the following adiabatic approximation, 



(g,a) 



aq 



27r crA; + aq — ie 



where 



(4.5) 



(4.6) 



Note that Eq. (|4.5p preserves the initial form of the spin- 
flip susceptibility given in Eq. (|3.16p , but with the initial 
gap Aq = (k^ — fcQ")/Ao is replaced by the flowing gap 
A; at scale I, and an overall reduction of the amplitude 
by the vertex correction jf. Indeed, using Eq. (j3.4ip we 
find Aj^o = Aq = {k^ — fco^)/Ao, so that for Z = we 
recover from Eq. (j4.5p the rescaled version of the initial 
condition (|3.16p . On the other hand, using the fact that 
the rf are fine tuned to reach a finite limit for Z — > oo, we 

oo. 
is 



see that A; ^ Af = e'(fc+ - fc")/Ao = e'A for / 



A justification for the adiabatic approximation 
given in the Appendix. 

To simplify the integrals, it is convenient to slightly 
modify the denominator in the expression for Tf{q, ie, a) 
in Eq. gH), 



aAr 



= ctA; + r? 



crA, 



(4.7) 



We have checked numerically from the solution of 
the resulting equations that this approximation is self- 
consistent by verifying the neglected term rf is indeed 
small. We thus arrive at the following approximation for 
the inhomogeneity rf{q,ie,a) that controls the flow of 
the fermionic self-energy. 



7i 



{2n) 

[F^*(Q,ze)]„oe"0+ 



1) 



i{€ + e) - a{q + q) - crA; 



(4.8) 



To be consistent the approximation (|4.7p we should also 
neglect the flowing anomalous dimension rj^ in the flow 
equation p.30p for rf , because an expansion of Eq. (|4.ip 
in powers of rf leads to a cancellation of the term rjfrf. 
The flow equation for rf then reduces to 



dirf = 



rr(o). 



(4.9) 



where 



7f[Fr(g,»e)]aa 
ie — aq — aAi 



Mr 



(4.10) 



Our general self-consistency equation p.4ip for the true 
Fermi point distance can then we written as an integral 
involving the flow of the couplings on the entire RG tra- 
jectory, 

/•OO 

A = Ao- / dZe"' Vcrrr(O) ■ (4.11) 



Anticipating that within our approximations rjf is inde- 
pendent of (7, we may write vif — rji. The flow equation 
for the spin-flip vertex then reduces to 



dai = -r]ai , (4.12) 
where from Eqs. (|3.2ip and (|4.8p we find 



dqd-e_^^^_^_^^2l^r^^^ (4.13) 



(2^)2 



[ie — aq — aAi]'^ 



We thus arrive at a closed system of flow equations for the 
two relevant couplings r^ and rf and the two marginal 
couplings gn,i and gd- We emphasize that our trunca- 
tion does not rely on a weak coupling expansion, which 
enables us to study a possible confinement transition. 

To give an explicit expression for the bosonic spin- 
flip propagator, we define dimensionless bare couplings 
(keeping in mind that we have chosen z/q = 1/vp), 

2j^oJaa = 27r5c,o , 2z/oJi;-a = 27rg„,o , (4.14) 
and the flowing couplings 

del = ifgcfi , gn,i = ifgnfi , (4.15) 

which according to Eq. (|4.12p satisfy the flow equations 

digc.i = -'2r]i9c,i , dign.i = -'2r]ign,i ■ (4.16) 
The rescaled spin-flip propagator can then be written as 



jf[friq,te)l 



1 - gc.o2^[nr (Q, +) + UfiQ, -)] + {gl, ~ .g2_o)(27r)2nr (Q, +)nf'^{Q, -) 
gc,i{ie + aq-aAi) + {gli-gli){aq~aAi) 



27r(ie — aq ~ crA; 



[ie-u}l{q,<7Ai) 



(4.17) 
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where 

uj^{q,x) = x{i - gc,i) ± '^x'^glj + - 5c,/)2 - glj 

= ^(1 - 9c.,i) ± ^x2(l - 5,,,)' + - - 9c,ir - aU ■ (4.18) 

Noting that W;^(0, x) = x{l — gc^iign,i)j we see that for small interaction strength both modes uj'^ and uj'^ are gapped. 
However, the gap of the mode i^/~(0, x) vanishes for ijc,; + = !> signaling a possible quantum phase transition to a 
confined state. In the present work we do not attempt to extend the RG beyond this point, but focus on the regime 
9c,i + gn,i < 1 were both modes are gapped. The integrals in Eqs. (|4.10p and (|4.13p can be carried out analytically 
using the residue theorem, with the result 

fr(0) = 2,5,,_ie(A, - l)5,,z + 9(1 - Az) 

and 

^ e (l - AQg^ 

^il-gc,r-gl,il^Af) 



9c,i 



(l-.9c,/)2-.9L(l 



A?) 



(4.19) 



1 - 9cj + ^ {I - gc,ir - gliil - Af) 



B. Self-consistent one- loop approximation 

The above system of coupled equations can only be 
solved numerically. However, if we neglect the flow of 
the coupling constants on the right-hand sides of these 
equations, we can obtain an approximate analytical solu- 
tion, which is useful to get a rough idea about the mech- 
anism responsible for the confinement transition. In this 
subsection we therefore set 



A/ 



9i,0 : 

A,* = Ae 



(4.21) 
(4.22) 



We expect that these approximations over-estimate the 
tendency towards confinement, because we know from 
Eq. (|4.16p that the flowing couplings gc,i and gn.i arc 
smaller than the bare ones. The second approximation 
(|4.22p is justifled provided the trajectory integral (|4.1ip 



is dominated by ? > 1 where A^ ^ 



Substituting 



Ae' on the right-hand side of Eq. (|4.1ip we obtain 



with 



Ao - A / dxl{x) 
J A 



i{x) = -e(x-i)% 



(4.23) 



-Q(l-x) 



a;y^(l-5c,o)2 -5^,o(l-a;2) 



(4.24) 



The a;-integration is elementary and we finally obtain 

Ao 



A 



1 + i?(A) 



(4.25) 



with 

R{A) = / dxl{x) 

J A 

•\/(l -3c,o)2 -5^,0 



In 



2ffc,o 
1 + 



9„,o 



A 1 



(l-9c,o)^-3i;.o 



(4.26) 



This expression diverges for g^^ + gn,o 1, correspond- 
ing to a confinement transition with A — )■ 0. Expanding 
i?(A) to second order in the couplings we obtain 

R{A) = -2g,.o + 2gl, ln(l/A) + 0{gl,) . (4.27) 

Keeping in mind that in Ref. [l3| we have neglected the 
chiral couplings and that here we have retained only in- 
terchain backscattering, Eq. (|4.27p is consistent with our 
previous weak-coupling result given in Eq. (jl.ip . From 
Eq. (|4.27p it is clear that there is a competition between 
the chiral part gc.o and the non-chiral part gnfi of the in- 
teraction. The chiral part gc,o leads to a repulsion of the 
Fermi points, while the non-chiral part gn,o generates an 
attraction and eventually triggers the confinement transi- 
tion, because for sufficiently small A the logarithm over- 
whelms the term linear in g^.a- 

To simplify the following analysis we shall from now 
on restrict ourselves to the special case gc,o — 0. This 
is sufficient to study the confinement transition, which is 
driven by the non-chiral part of the interaction. Setting 
<?o = 9nfl the confinement transition occurs within the 
approximations in this section at go = 1- A numerical 
solution of Eqs. (|4.25p and (|4.26p for the true A as func- 
tion of go is shown in Fig. [5] The confinement transition 



11 



A/Ao 




FIG. 6: Numerical solution of Eq. (I4.25|) for gcfi = as func- 
tion of go = Qnfi for different values of Ao — {k^ — fc[5')/Ao. 



for 50 ^ 1 is clearly visible. In fact, the behavior of A 
for go ^ 1 can be obtained analyticaUy. In this case 
A <C V 1 — ffo ' so ti^sA, we may approximate 



R{A) 



:ln 



A 



(4.28) 



The self-consistency condition (|4.25p for A then reduces 
to 



A; 



Ao 



2 In 



2V1^ 
A 



(4.29) 



For A <C ■\/l — (?o the second factor in the square braces 
of Eq. (|4.29p is small compared with unity, so that it is 
consistent to take the limit go I in this expression. If 
we identify A with the order parameter of the confine- 
ment transition transition (with A corresponding 
to the deconfined phase), then Eq. (|4.29[) predicts mean- 
field behavior with logarithmic corrections. 

Our simple one-loop approximation thus predicts that 
for go —^ 1 there is a confinement transition where the 
true Fermi point distance A collapses, corresponding to 
vanishing effective interchain hoping t*^ ^ 0. In pseudo- 
spin language, the quantum critical point go = 1 cor- 
responds to a vanishing magnetization in z-direction, in 
spite of the fact that there is a uniform magnetic field. 
For go > 1 our one-loop approximation suggests that 
there is long-range ferromagnetic order in xy-direction. 
However, in one dimension we do not expect true long- 
range order, so that fluctuations beyond the one- loop ap- 
proximation should be important. It is therefore impor- 
tant to go beyond this approximation, which we shall do 
in the following subsection. 



C. Including the renormalization of the effective 
interaction 

We now improve the above calculation by taking the 
flow of the effective interaction into account. For sim- 
plicity, we focus again on the special case without chiral 
interactions, so that we need to keep track only of the 
flowing non-chiral interaction gn^i = gi- Furthermore, 
in the absence of chiral couplings rf — ari. The self- 
consistency equation (|4.11[) for the Fermi point distance 
then reduces to 



Ao 



die' 



29(1 - A;)A,.g2 



where 



Ai^At- 2ri = Ae' - 2n 



and the flow of r; and gi is determined by 
diri = ri +A{gi,Ai) , 



digi = B{guAi) 



with 



and 



A{guAi)^ 


Q{l-Ai)Aig^ 


^1-5^(1-^?)' 


B{gi,Ai) - -2r]igi 

-26(1 - AOfff 




l + \/l-5f(l-A2) 



(4.30) 

(4.31) 

(4.32) 
(4.33) 

(4.34) 



(4.35) 



The initial value has to be fine tuned such that 
lim;_»oo ri remains finite. This leads to the self- 
consistency equation (j4.30p for the true Fermi point 
distance. We emphasize again that our approxima- 
tion scheme is not based on a weak coupling expan- 
sion, so that the /3-function given in Eq. (|4.35p is non- 
perturbative in the coupling gi. Instead of Eq. (|4.26[) we 
now obtain for the dimensionless renormalization factor 
R{K) defined in Eq. g^S), 



i?(A) = 



2 
A 



dle-'A{guAi) 



2 / ' dl{l - 2e-V,/A)^£L:^M= . (4.36) 

Note that formally the confinement transition manifests 
itself via a divergence of the function i?(A) for A — > 0. 
However, as shown in Fig. [7l the renormalization factor 
remains now finite so that also the self-consistent A is 
finite for go = 1, as shown in Fig. [51 We conclude 
that the confinement transition obtained in the previous 
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0.2 0.4 0.6 0.8 1 



go 

FIG. 7: Numerical solution of the renormalization factor 
R{A) defined in Eq. (|4.36p as a function of the bare coupling 
go for different values of Ao. 



A/A, 




I ' ' ' ' ' ' ' ' ' 1 

0.2 0.4 0.6 0.8 1 

go 



FIG. 8: Numerical solution of the true Fermi point distance 
A = (fc^ — fc~)/Ao as a function of the bare coupling go for 
different values of the bare distance Ao = {k^ — kQ)/Ao. 



subsection is an artefact of the approximations (|4.2ip and 

Km . 

Let us take a closer look at the point go — ^ where 
the one-loop approximation predicts a confinement tran- 
sition. The RG flow of the couplings gi and r; as well 
as the flowing anomalous dimension -qi for this case is 
shown in Fig. [SI One clearly sees that the initially large 
flowing anomalous dimension drives the running coupling 
gi towards smaller values; eventually gi approaches a 
finite limit for large I. Furthermore, the running cou- 
pling ri approaches its asymptotic value for large I non- 
monotonously. The true Fermi point distance A as a 
function of the bare one for go = 1 is shown in Fig.lTOl We 
see that large interchain backscattering strongly reduces 
the Fermi point distance, although the Fermi surface 
never collapses, in agreement with scenario suggested by 
the weak coupling analysis^^. 





— 9i / 















0.2 0.4 0.6 0.8 1 

l/ln(1/A) 



FIG. 9: RG flow of the couplings gi, ri and the anomalous 
dimension i]i as a function of the logarithmic flow parameter 
I for Ao = 0.1 and go = 1. 
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0.2 0.4 0.6 0.8 1 

Ao 

FIG. 10: Numerical solution of the true Fermi point distance 
A for go — 1 as a function of the bare distance Aq. 

V. CONCLUSIONS 

In this work we have used a functional RG approach to 
calculate self-consistently the true distance A — — 
between the Fermi points of the bonding and the an- 
tibonding band in a system consisting of two chains of 
spinless fermions connected by weak interchain hopping 
tj_. Using the insight from our earlier weak coupling 
analysisi^ that the renormalization of the Fermi surface 
is essentially determined by interchain backscattering, we 
have treated this scattering process non-perturbatively 
by representing it in terms of a collective bosonic field x- 
In pseudospin language, where t± = h corresponds to a 
uniform magnetic field in z-direction and the interchain 
backscattering interaction corresponds to a ferromag- 
netic xy-interaction, the field x can be viewed as a fluc- 
tuating transverse magnetic fleld, which competes with 
the uniform field h in z-direction. A self-consistent one- 
loop approximation predicts that for sufficiently strong 
interchain backscattering there is indeed a quantum crit- 
ical point were the renormalized Fermi point distance 
A (X t*j_ vanishes. However, a more accurate calculation 
taking vertex corrections and wave-function renormaliza- 
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tions into account shows that the renormahzed A (x t*j_ 
remains finite. This is in agreement with the expectation 
that a ferromagnetic ccy-interaction in a one-dimensional 
itinerant electron gas cannot give rise to long-range fer- 
romagnetic order. Previous studies of the spinless two- 
chain proble m^^'^^ came to the conclusion that the sys- 
tem exhibits a confinement transition if the anomalous 
dimension rjQ of the Luttinger liquid for ij^ = is unity. 
The important difference between these earlier works and 
our calculations is that we have completely neglected pair 
tunneling. In a subsequent article® we shall show how 
the inclusion of this process stabilizes again the flow of 
the interchain backward scattering and enhances the ten- 
dency towards confinement. In the same article we shall 
also show that our approach can be generalized to study 
the more interesting and physically more relevant con- 
finement problem in an infinite array of weakly coupled 
metallic chains. In this case the Fermi surface consists 
of two disconnected sheets, which self-consistently de- 
velop completely fiat sectors at the confinement transi- 
tion. We have preliminary evidence that in this case there 
exists a confined phase where the renormahzed interchain 
hopping vanishes. The essential scattering process driv- 
ing this transition is the non-chiral part of the density- 
density interaction which transfers momentum within a 
given sheet of the Fermi surface. 

Finally we point out that this work describes also some 
technical progress: we have been able to find a sensible 
extrapolation of the weak coupling functional RG ap- 
proach to the strong coupling regime. Our truncation 
strategy of the formally exact hierarchy of functional RG 
equations relies on the classification of the vertices ac- 
cording to their relevance. An obvious disadvantage of 
this approach is that we cannot give reliable error esti- 
mates, which is a common feature of most truncations 
of the coupled functional RG flow equations for the ver- 
tex functions. Note, however, that a similar truncation 
of the functional RG equations for the interacting Bose 
gas gave quite accurate results for the shift in the critical 
temperature^. In the present problem we know a pri- 
ori from the weak coupling analysis that the physics is 
dominated by a single scattering channel, the inter-chain 
backscattering. However, an extension of our approach 
to problems where several scattering channels compete 
seems to be possible. 
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FIG. 11: Additional diagrams contributing to the flow of the 
spin-flip susceptibility in a cutofl scheme with a fermionic 
band-width cutofl. The solid arrow with slash is the fermionic 
single-scale propagator. 



APPENDIX: JUSTIFICATION OF THE 
ADIABATIC APPROXIMATION 



We give here a justification for the adiabatic approxi- 
mation for the rescaled spin-fiip susceptibility given in 
Eq. (|4.5p . Let us therefore use a more general two- 
cutoff procedure where we impose a band-width cutoff 
A^ — A^e""^ on the fermionic propagator in addition to 
the bosonic momentum transfer cutoff A; = Aqc^K All 
vertices and coupling constants then depend on both log- 
arithmic fiow parameters I and r. Instead of Eq. (|3.35p 
the rescaled fermionic propagator can then be approxi- 
mated by 



e(l< |g| <eO 



and the corresponding single-scale propagator is 

Si\q\-1) 



It - avi .^q - ri ,^ 



(A.l) 



(A.2) 



We recover the vertices of the momentum transfer cutoff 
scheme by taking first the limit r oo. Of course, the 
result of the RG should be independent of how we reach 
a certain point in two-dimensional cutoff space spanned 
by A; and Af . Suppose we first fix Af = A^ and 
perform the reduction of the momentum transfer cut- 
off Aq — > A/. As a second step we reduce the fermionic 
cutoff A^ 0. On the right-hand side of the fiow equa- 
tion for the spin-fiip susceptibility there are then two ad- 
ditional diagrams^ involving the fermionic single-scale 
propagator, which are shown in Fig. 1111 For a fixed scale 
I the corresponding fiow equation is (for simplicity we set 
a = I and omit the chirality label) 



-[qd^ + ed,]iiiUQ) + ni:?(Q),(A.3) 
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with 



~Gl,iQ)GlAQ + Q + aA:) 



(A.4) 



where A* = A/Af . The flow equation for the spin-flip 
vertex is now of the form 



dr^l,r^e{X + l~T)Cil,T) , 



(A.5) 



where C{1,t) is some function of the flow parameters and 
of the running couphng constants, and A = ln(A^/Ao). 
For simphcity we choose Aq — Aq , so that A = 0. The 
0-function is due to the fact that the internal foop mo- 
menta are restricted by the momentum transfer cutoff 
Aq. Obviously, dr^i^r = for r > /, so that 7; = 7/,r>i is 
independent of r. Similarly, the flow equation for rf^ is 
of the form 



drrrr=rl, + eil-r)A{l,T) , 



(A.6) 



with some other function j4(/,r). This implies rf^ — 
e^^^rf^^i for r > Z, so that the flowing Fermi point dis- 
tance A;^T-7 defined analogous to Eq. (|4.6p via 



A,. 



Af 



1,tI 



(A.7) 



is for r > ^ of the form 



A _ ^ 



(A., 



' l.T=l 



and A/ = A/A,-(r^ 



-I- 



where we have defined rf = 
Tj"), see Eq. dm]). For t — > 00 the solution of Eq. (|A.3|) 
can therefore be written as 



nr(Q) 



dT tlfriQe^^'iAie^-') 



(A.9) 



Using the fact that in the integral of the last term we may 
pull out a factor of and approximating the fermionic 
Green functions by Eqs. (jA.ip and (jA.2[) with vf^ and 
n,T set equal to zero, we recover from the last term the 
adiabatic approximation (|4.5p . Actually, to calculate the 
flow of the fermionic self energy, we only need the bosonic 
Green hmction at 5 = ±1. For A; < 1 the non- adiabatic 
contribution to the polarization then vanishes, leaving 
us just with the adiabatic part. For larger A; there is a 
crossover to the more general expression. However, the 
leading contribution to the flow of A; itself stems from 
the region where the adiabatic approximation is valid. 

Finally we note that the adiabatic approximation (14. 5p 
is also consistent with the flow equation for the spin-flip 
susceptibility in the momentum transfer cutoff scheme 
shown in Fig. [3J taking the derivative of the right-hand 
side of Eq. (|4.5p with respect to the flow parameter and 
inserting for the derivative of the fermionic self-energy 
its flow equation (|3.5p . we find that the right-hand side 
can be written in terms of a symmetrized closed fermion 
loop with four bosonic legs and renormalized propaga- 
tors. This corresponds to the adiabatic approximation 
for the vertex F^'^-* which according to Fig. [3] drives the 
flow of the spin-flip susceptibility. 
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